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In this work, we consider time-dependent dark energy star models, with an evolving parameter 
^ • uj crossing the phantom divide, ui — — 1. Once in the phantom regime, the null energy condition 

is violated, which physically implies that the negative radial pressure exceeds the energy density. 
•"T' , Therefore, an enormous negative pressure in the center may, in principle, imply a topology change, 

consequently opening up a tunnel and converting the dark energy star into a wormhole. The criteria 
\Q ' for this topology change are discussed and, in particular, we consider a Casimir energy approach 

involving quasi-local energy difference calculations that may reflect or measure the occurrence of 
a topology change. We denote these exotic geometries consisting of dark energy stars (in the 
phantom regime) and phantom wormholes as phantom stars. The final product of this topological 
change, namely, phantom wormholes, have far-reaching physical and cosmological implications, as in 
addition to being used for interstellar shortcuts, an absurdly advanced civilization may manipulate 
these geometries to induce closed timelike curves, consequently violating causality. 
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Q ■ I. INTRODUCTION 
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Recent high-precision observational data have confirmed that the Universe is undergoing a phase of accelerated 
expansion Several candidates, responsible for this expansion, have been proposed in the literature, in particular, 
dark energy models (see Ref. [2| for a review) and modified gravity (e.g., see Refs. Q for recent reviews). In 
particular, the former models are fundamental candidates, in which a simple way to parameterize the dark energy 
is by an equation of state of the form ui = p/p, where p is the spatially homogeneous pressure and p is the dark 
energy density. A value of ui < —1/3 is required for cosmic expansion, and u> = — 1 corresponds to a cosmological 
constant. A specific exotic form of dark energy denoted phantom energy, with ui < — 1, has also been proposed 0], 
and possesses peculiar properties, such as the violation of the null energy condition (NEC) and the energy density 
increases to infinity in a finite time at which point the size of the Universe blows up in a finite time, which is 
known as the Big Rip. In this context, the violation of the NEC presents us with a natural scenario for the existence 
of traversable wormholes, and indeed it has been shown that these exotic geometries can be supported by phantom 
energy [H, [||. It is also interesting to note that recent fits to supernovae, CMB and weak gravitational lensing data 
probably favor an evolving equation of state, with the parameter crossing the phantom divide ui = — 1 

Despite the fact that the dark energy equation of state represents a spatially homogeneous cosmic fluid and is 
assumed not to cluster, it is possible that inhomogeneities may arise due to gravitational instabilities. More precisely, 
although the equation of state leading to the acceleration of the Universe on large scales is an average equation 
of state corresponding to a background fluid, it is possible that dark energy condensates may possibly originate 
from density fluctuations in the cosmological background, resulting in the nucleation through the respective density 
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perturbations. Despite the fact that once in the dark energy regime the material system becomes gravitationally 
repulsive, we may consider the possibility of the formation of a matter system that originally obeys all the energy 
conditions. Cosmological observations do not rule out, and in some studies favor, an evolving equation of state for 
the dark energy. It is therefore quite possible that what we know as dark energy today has evolved from a more 
benign fluid. An over-density of this fluid could in principle commence a collapse into a star which. Such a model is 
presented in section UlI CI We also point out that even in the case of a dark-energy fluid, there is no definite resolution 
to the debate of clustering scales. This is mainly due to non-linearity, especially in the vein of dark energy interacting 
with ordinary fluids. It may also be possible to glean some information on the cosmological dark matter by studying 
certain properties of such gravitational condensates. (See [8[ and references therein for comments on these issues.) In 
this context, a number of inhomogeneous solutions have been the object of analysis, such as the phantom wormholcs 

[f| mentioned above, dark energy stars @, and other structures such as condensates supported by the generalized 
Chaplygin gas (Io| which possibly arise from density fluctuations in the generalized Chaplygin gas background, and 
condensed structures supported by the van der Waals equation of state [llj • In a recent paper [12} , it was also shown 
that the 4D Einstein-Klein- Gordon equations with a phantom scalar field possess non-singular, spherically symmetry 
solutions, although a stability analysis on these solutions indicates they are unstable. 

The dark energy star models are also a generalization of a new emerging picture for an alternative final state of 
gravitational collapse, namely, the gravastar (gravitational vacuum star) models. The latter proposed by Mazur and 
Mottola [l3j . has an effective phase transition at/near where the event horizon is expected to form, and the interior 
is replaced by a de Sitter condensate. The latter is then matched to a thick layer, with an equation of state given by 
p = p, which is in turn matched to an exterior Schwarzschild solution. The issue of gravastars has been extensively 
analyzed in the literature, and we refer the reader to Refs. [H], EH]. The generalization of the gravastar picture is 
considered by matching an interior solution governed by the dark energy equation of state, lo = p/p < —1/3, to an 
exterior Schwarzschild vacuum solution at a junction interface Q. The dynamical stability of the transition layer 
was also explored, and it was found that large stability regions exist that are sufficiently close to where the event 
horizon is expected to form, so that it was argued that it would be difficult to distinguish the exterior geometry of 
the dark energy stars from an astrophysical black hole. Thus, these alternative models do not possess a singularity at 
the origin and have no event horizon, as its rigid surface is located at a radius slightly greater than the Schwarzschild 
radius. This restriction arises from the observed lack of energy emission due to surface collisions of infalling material 
in suspected black hole systems. In fact, although evidence for the existence of black holes is very convincing, a certain 
amount of scepticism regarding the physical reality of event horizons is still encountered, and it has been argued that 
despite the fact that observational data do indeed provide strong arguments in favor of event horizons, they cannot 
fundamentally prove their existence [l6j |. 

As mentioned above, recent fits to observational data probably favor an evolving equation of state, with the dark 
energy parameter crossing the phantom divide u> = — 1 [7|. Motivated by this fact, in a rather speculative scenario 
one may theoretically consider the existence of a dark energy star, with an evolving parameter starting out in the 
range — 1 < uj < —1/3, and crossing the phantom divide, u) = — 1. Once in the phantom regime, the null energy 
condition is violated, which physically implies that the negative radial pressure exceeds the energy density. Therefore, 
an enormous negative pressure in the center may, in principle, im ply a topology change, consequently opening up 
a tunnel, and converting the dark energy star into a wormhole Il7|. One may assume that the topology change 
may occur at approximately the Planck length scales, and once created may be self-sustained as shown in Ref. 
[18} . In fact, the change in topology is an extremely subtle issue, as in general relativity these changes probably 
entail spacetimc singularities. However, at the Planck length scales quantum gravity effects dominate and spacetime 
undergoes a deep and rapid transformation in its structure, probably producing a multiply-connected quantum foam 
structure 0, [2(| ■ It was suggested in Ref. [TtJ that one could imagine an absurdly advanced civilization 21 1 pulling 
a wormhole from this submicroscopic spacetime quantum foam and enlarging it to macroscopic dimensions. However, 
in a more plausible scenario, the possibility that inflation might provide a natural mechanism for the enlargement of 
such wormholes to macroscopic size was explored [22[. In this work, we outline the theoretical difficulties associated 
to the change in topology and present a method based on the Casimir energy approach. Although it is still unsure 
if this method produces a topology change, it is extremely useful as the quasi-local energy difference calculation may 
reflect or measure the occurrence of a change in topology. Other concepts of to pol ogy changing spacetimes have been 
studied, for instance: Using semi-classical and Morse-index methods in Refs. [23]; higher order back-reaction terms 
due to fluctuations of gauge fields in the vicinity of a black hole may result in the formation of a wormhole-like object 
[24} ; and more recently an approach based on a Ricci flow may result in quantum wormholes [2f| . 

Once the topology change has occurred, with the respective opening of a tunnel, then the dark energy star has 
been converted into wormhole supported by phantom energy. As mentioned above, it has recently been shown that 
traversable wormholes may, in principle, be supported by phantom energy 0, which apart from being used as 
interstellar shortcuts, may induce closed timelike curves with the associated causality violations [26|, [27| • Particularly 
interesting solutions were found @ , and by using the "volume integral quantifier" , it was found that these wormhole 
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geometries are, in principle, sustained by arbitrarily small amounts of averaged null energy condition (ANEC) violating 
phantom energy. A complementary approach was traced out in [5| , by considering specific choices for the distribution 
of the energy density threading the wormhole. Recently, AD static wormhole solutions supported by two interacting 
phantom fields were found as well 28]. Despite the fact that traversable wormholes violate the NEC in general 
relativity (see Ref. [29( for a recent review), it has been shown that the stress e nerg y tensor profile may satisfy 
the energy conditions in the throat neighborhood in dynamic wormholes (see Ref. [30] and references therein) and 
in certain alternative theories to general relativity [311 ]. Perhaps not so appealing, one could denote these exotic 
geometries consisting of dark energy stars (in the phantom regime) and phantom wormholes as phantom stars. We 
would like to state our agnostic position relatively to the existence of dark energy stars and phantom wormholes, or 
for that matter of phantom stars. However, it is important to understand their general properties and characteristics, 
and we emphasize that the presence of a dark energy fluid permeating the universe makes the study of dark energy 
condensates a physically relevant endeavor. 

This paper is organized in the following manner: In section UH we briefly review static dark energy stars, followed 
by a deduction of general solutions of time-dependent spacetimes. In Section ITTTl specific time-dependent dark energy 
solutions are outlined, in particular, we present the specific cases of a constant energy density, the Tolman-Matese- 
Whitman mass function solution, and a class of models with a non-zero energy flux term, which form from gravitational 
collapse. In Section IIV1 we describe the theoretical difficulties associated with changes in topology and present in 
some detail specific methods used in the literature, namely, a Casimir energy approach involving quasi-local energy 
difference calculations that may reflect or measure the occurrence of a topology change. We also briefly review the 
Morse index analysis. In section [V] we conclude. 



II. TIME-DEPENDENT DARK ENERGY STARS 



A. Static space-time 



In this section, we provide a brief outline of the mathematical models of static and spherically symmetric dark 
energy stars considered in Ref. Consider the following time-independent line element, in curvature coordinates, 
representing a dark energy star 

rlr 2 

ds 2 = ~e 2a ^ dt 2 + - , w + r 2 (d9 2 + sin 2 B dcp 2 ) , (1) 

1 — 2m(r)/r 

where a(r) and m(r) are arbitrary functions of the radial coordinate, r. The function m(r) can be interpreted as the 
quasi-local mass, and is denoted as the mass function Q. The factor a(r) is the "gravity profile" and is related to 
the locally measured acceleration due to gravity, through the following relationship: A = \J\ — 2m(r)/ra'(r) [1, [Il|, 
where the prime denotes a derivative with respect to the radial coordinate r. The convention used is that a'(r) is 
positive for an inwardly gravitational attraction, and negative for an outward gravitational repulsion. 
The Einstein field equations are given by @ 

m ' = 47rr 2 / 9, (2) 
, m + 4irr 3 p r 



r(r — 2m) 



(3) 



(p + p r ){m + 47rrV) 2 

Pr = 7 + -{Pt- Pr) , (4) 

r(r — 2m) r 

where p(r) is the energy density, p r (r) is the radial pressure, and pt(r) is the tangential pressure orthogonal to p r . 
Note that Eq. (01 corresponds to the anisotropic pressure Tolman-Oppenheimer-Volkoff (TOV) equation. 

An additional constraint is placed on the system of equations by considering the dark energy equation of state, 
p r (r) — u>p(r), and taking into account Eqs. @ and ([3]), we have the following relationship 

, m + Lurm' 

a ( r ) — — ; ^ — r • (5) 

w r (r - 2m) w 

There is, however, a subtle point that needs to be emphasized 0, [f|. The notion of dark energy is that of a spatially 
homogeneous cosmic fluid. Nevertheless, it can be extended to inhomogeneous spherically symmetric spacetimes, by 
regarding that the pressure in the equation of state p = up is a radial pressure, and that the transverse pressure 
may be obtained from Eq. In addition to this, and as mentioned in the Introduction, despite the fact that 
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the dark energy equation of state represents a spatially homogeneous cosmic fluid and is assumed not to cluster, 
inhomogcncities may arise due to gravitational instabilities. Thus, the dark energy star geometries considered here 
may possibly originate from density fluctuations in the cosmological background, resulting in the nucleation through 
the respective density perturbations fill ] - 

In Ref. Q , specific solutions were found by considering that the energy density is positive and finite at all points in 
the interior of the dark energy star. In particular, several relativistic dark energy stellar configurations were analyzed 
by imposing specific choices for the mass function m(r), and through Eq. ([5]), a(r) was determined, consequently 
providing explicit expressions for the stress-energy tensor components. This interior solution was further matched to 
an exterior Schwarzschild vacuum solution given by 



ds 2 = - 1 ) dt 2 + I dr 2 + r 2 (d9 2 + sin 2 d(f> 2 ) , (6) 



2M\ , / 2M\~ 1 
-) dt+ { 1 -—) ' 

at a junction interface a. The Schwarzschild spacetime possesses an event horizon at = 2M , so that to avoid the 
latter, the junction radius lies outside 2M, i.e., a > 2M. 
The surface stresses on the thin shell are given by 

I oZ \ 

(7) 
(8) 

where a and V are the surface energy density and the tangential surface pressure [Hil HI, respectively The overdot 
denotes a derivative with respect to t, which is the proper time on the junction interface, and the prime here denotes 
a derivative with respect to the junction surface radius a. 

The dynamical stability of the transition layer a of these dark energy stars to linearized spherically symmetric radial 
perturbations about static equilibrium solutions was also explored. It was found that large stability regions exist that 
are sufficiently close to where the event horizon is expected to form, so that it would be difficult to distinguish the 
exterior geometry of the dark energy stars, analyzed in 0], from an astrophysical black hole. 




B. Time-dependent spacetime 



In this section, we generalize the above static dark energy star models to time-dependent geometries. This is mainly 
motivated by the fact that recent fits to supernovae, CMB and weak gravitational lensing dataprobably favor an 
evolving equation of state, with the dark energy parameter crossing the phantom divide u> = — 1 [7[. 

In the following, we consider a time-dependent and spherically symmetric metric given by 



ds 1 



Note that one may also define the function f3(r, t) as 

/?(r,t) = ; 



■In 



dr 1 



1 - 



(de 2 



2m(r, t) 



(9) 



(10) 



where the 'mass function' m(r, t) is now time-dependent. 

The Einstein field equation provides the following nonzero components: 



P(r,t) 
p r (r, t) 

f(r,t) 

Pt(r, t) 



-2/3 



8irr 2 



-2;1 



{2(3' r - 
8*r>-< 2aV - 



1) 



1), 



47rr 



8tt' 



-2/3 



a + a 



12 



a' 13' + -(a'- 13') 
r 



+ e~ 2a (a(3-P- P 2 ) \ , 



(11) 
(12) 

(13) 
(14) 
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where the prime denotes a partial derivative with respect to the radial coordinate r, and the overdot a partial derivative 
with respect to the time coordinate t. Note the presence of an energy flux term in the radial direction, T r = ±/(r, t), 
which depends on (3. 

An important issue in the time-dependent dark energy stars with the parameter crossing the phantom divide are 
the energy conditions, in particular, the null energy condition (NEC). The NEC is defined as T^k^k" > 0, where 
k^ is any null vector, and consequently provides p + p r ± 2/ > 0. The latter definition, taking into account the field 
equations (fTTj) - ([T3l) , is given by 



p + Pr ±2f = 



1 



e" 2/ V + /3') ± 2/3e-( Q+ « 



> 0. 



In this context, one may consider a generalization of the equation of state p r — u>p, given by 

Pr {r,t) =w(r,t) [p(r,t)±2f(r,t)} , 



(15) 



(16) 



p r (r,t) T2f(r,t) = uj(r,t)p(r,t) 



(17) 



where f(r,t) — ^T 1 ^. is the energy flux term, as noted above. However, one can come up with an interesting class of 
solutions considering the following equation of state: 



p r {r,t) = uj(r,t)p(r,t) 



(18) 



Throughout this work, we essentially use the equation of state given by Eq. (118)) . 

Taking into account the field equations (lll|l and (fT2|) . then Eq. (fT8|) provides the following relationship 



w(r, t)p'(r, t) = a'(r, t) + —[l + w(r, t)} 

2r 



1 



,20{_r,t) 



A similar analysis was carried out in Ref. [34|. in the context of time-dependent wormholcs. 

Equation (|19| may be formally solved in terms of a(r,t), and provides the following general solution 



(19) 



a(r,t)= / -\u{r,t)0'(f,t)r--[l + Lj(f,i)] 



1 - e 



2/3(f,t) 



dr . 



Thus, in principle, if uj(r, t) and (3{r,t) are known, then a(r, t) may be obtained from Eq. (|20p . 
As an alternative, Eq. (fl9|) can be formally integrated for f3(r,t) to yield the general solution 



/?(r,t) = -~ln 



-2i?(t) 



e r (^)(l +uj(f, t)) 



w(r, t)r 

where F(t) is an integration function, and the factor T(r, t) is defined as 

"2ra'(f,i) + f+l 



r(r,t) 



roj(r, i) 



(20) 



(21) 



(22) 



A particularly simple and interesting toy model is the specific case of a purely time dependent parameter to = u>(t), 
so that the general solution (|2"Tj) takes the form 



l/ 2a(r,t) 
2 | w(i) 



-In 



uj(t) 



(23) 



In the next section, we analyze specific solutions, namely, that of a constant energy density, the Tolman-Matese- 
Whitman mass function, which were extensively explored in Ref. Q, and a collapsing model with a nonzero energy 
flux term. 
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III. SPECIFIC TIME-DEPENDENT DARK ENERGY STAR SOLUTIONS 

A. Constant energy density 

Consider the specific case of a constant energy density, p(r, t) — po, so that Eq. (fTTj) provides the solution 

2 , Fi(ty 



1-Ar* 



(24) 



with A = 87T/9o/3 and F±(t) is a function of integration. 
Substituting into Eq. (|2"0|) , one arrives at 



a(r, t) 



3Auj(r,t)r 3 + Ar 3 - 3F x (t) 
2r [3r- Ar 3 +3F 1 {t)] 



(25) 



where F 2 {t) is another function of integration. One may further simplify the analysis by considering that F\(t) = 0, 
which is physically justified by the imposition of a finite mass function at the origin r = for all values of the time 
coordinate t. Note that considering F±(t) = implies that the mass function is not time-dependent, and the flux term 
f(r,t) is zero, as $ = 0. 

For instance, consider the specific example of a separation of variables of the parameter w(r, t) given by 



Choosing the following functions: 



u(r,t) = u>i(t)oj 2 {r) 



LJi(t) = ujq + luq tanh [a(t — to)] , 
A 2 



w 2 (r) 



1 + {r/R) 2 ' 



(26) 

(27) 
(28) 



where ujq, ujq, a, t$, A 2 , and R are constants. The factor 1/er may be interpreted as the "relaxation time", describing 
the width of the time-dependence. 

See Fig. [T]for a qualitative description of w(r, t) given by Eqs. (|26"|) - (j28|) . The left plot represents the behavior of 
lo = Ld{0,t) at the center r = 0. We have considered the following values: R = 1, 6 = 1, Wq = 3/4, Q — 2/3 and 
to = 8. 
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FIG. 1: Left plot: For this case we considered the behavior of u — uj(0,t) at the center r = 0; for the dashed curve we have 
a — 0.3 and for the solid curve a — 0.1. Right plot: For this case we considered the uj — u)(t, r) dependence with a = 0.3. For 
both cases, we have assumed the following numerical values: R = 1, 5 = 1, ujq = 3/4, <Do = 2/3 and to = 8. 



Substituting the functions (|26p - (j28|) into Eq. (|25|) . yields the following solution 



a(r,t) 



In. 



(1 - Ar 2 )^ 1 ^' 1 
[l + (r/R) 2 ] AuJl{t) 



(29) 
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where, for notational simplicity, the constant A is defined as 

A 3AR2X2 (30) 

The function F%(t) can be absorbed through a redefinition of the time coordinate as before, so that without a significant 
loss of generality one may impose the condition F 2 (t) = 0. 
The pressure profile is given by 

3X 2 R 2 Au; 1 (t) 

Pr ~ ^{R 2 +r 2 ) > W 

_ JL4 R 2 Ar 2 (R 2 + 2r 2 ) + Ar 6 + u; 1 (t)X 2 R 2 [3Ar 2 u Jl (t)X 2 R 2 - A{Ar A + R 2 )] 
Pt ~32^ (1- Ar 2 )(R 2 + r 2 ) 2 ' ( ' 

with p r = pi at the center, r = 0. 

The analysis simplifies by considering a purely time-dependent parameter, i.e., lj = w(t). Thus, Eq. (|25p takes the 
form 

a{r,t) = -~[l + 3u/(t)]hi(l- Ar 2 ) + F 2 (t) . (33) 

The factor can be absorbed into a redefinition of the time coordinate, so that without a significant loss of 

generality, one can assume F 2 (t) = 0. 

The pressure profile is given by the following relationships: 

3Aj(t) 

ZA[^(t)+Ar 2 +ZAr 2 u 2 (t)] 

Pt = 32.(1 - Ar 2 ) ' (34) 

Note that p r — p t at the center, r = 0, as expected. 



B. Tolman-Matese- Whit man mass function 



An interesting example is the Tolman-Matese- Whitman mass function considered in Ref. [9(. As in the example 
outlined above, we impose that $ = 0, so that the flux term f(r, t) is zero. Thus, consider the following choice for the 
time-independent mass function, given by 



(35) 



where bo is a non-negative constant The latter may be determined from the regularity conditions and the finite 
character of the energy density at the origin r = 0, and is given by bo = 87rp c /3, where p c is the energy density at 
r = 0. 

Now, consider the radial and temporal dependent case of u> = u>{r,t) given by the functions (f26]) - (|28p . Substituting 
these functions and Eq. ([35]) into Eq. (|20|) . yields the following solution 



a(r,t) = ^ln{(l + b r 2 f {t) (l + 2bor 2 ) 



2\T(t) 



2(0' 



+ F 2 (t), 



(36) 



where the F 2 (t) is a function of integration which may be reabsorbed in a redefinition of the time coordinate, so that 
without a loss of generality we impose F 2 (t) = 0, as before. For notational simplicity, we have considered the following 
definitions 



E(t) 
T(t) 
S(t) 



1 + b R 2 (2b R 2 - 3) - b Q \ 2 R 2 uj 1 (t)(l - 2b Q R 



2(l-2R 2 b )(l-R 2 bo) 



2b X 2 R 2 
1 - 2R 2 b ' 

bo\ 2 R 2 m 1 (t)(2boR 2 -3) 
2(1-2E 2 6 )(1 -^6 ) ' 



(37) 
(38) 
(39) 
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respectively. 

The stress-energy tensor components are given by 
b {3 + 2b r 2 ) 



8ir(l + 2b r 2 ) 2 ' 
b R 2 X 2 {3 + 26 r 2 )wi(i) 



p(r) 

Pr ^~ 8tt(R 2 +r 2 )(l + 2b Q r 2 ) 2 ' 
Pt (r, t) = b Ui? 4 6gr 6 (l - wi(t)A 2 ) 2 + 8i? 2 &gr 8 (l + ^i(t)A 2 ) + 4i? 4 & 2 r 4 (2 - 3u>i(i)A 2 + 3^!(i)A 4 ) 
+ 4i? 2 6 2 r 6 (4 + 9cJi(i)A 2 ) + R%r 2 {3 - 16wi(i)A 2 + 9w 2 (t)A 4 ) + 2i? 2 6 r 4 (3 + 14 Wl (t)A 2 ) 
+ 6 r 6 (3 + 8fe r 2 + 46 2 r 4 ) - 12wi(i)i? 4 A 2 ] / [32tt(1 + b r 2 )(l + 2b r 2 ) 3 (R 2 + r 2 ) 2 ^ 



(40) 
(41) 



(42) 



Note that p r = pt at the center, r = 0. 

For simplicity, considering a purely time-dependent parameter to = u>(t), and substituting (|35[) into Eq. (|20[) . 
provides the following solution 



(1 + 6 r 2 ) 2 (1 + 26 r 2 ) 



f(t) 



(43) 



where the F(t) is a function of integration which, as before, may be absorbed into a redefinition of the time coordinate, 
so that one may consider F(t) = without a significant loss of generality. 
The stress-energy tensor components are given by 



p r = ui{t)p 



6 (3 + 2b r 2 )w{t) 



Pt 



8tt(1 + 2& r 2 ) 2 ' 

fr {46 2 r 4 (l + cj)[3 + 6 r 2 (l + u)] + 6 r 2 [3 + uj{9uj + 16)] + 8b 2 r 4 + 12lj} 
327r(l + 6 r 2 )(l + 26 r 2 ) 3 



(44) 
(45) 



with p r = p t at the center, r = 0. 



C. A class of models with a non-zero flux term 



In this section, we construct a set of models with a non-zero energy flux term, where at early times possesses a small 
inhomogeneity in the region near r = which grows due to gravitational collapse. For the specific case considered in 
this section, we assume for simplicity that the parameter, which eventually crosses the phantom divide in the central 
region, is purely time-dependent, i.e., to — Lu(t), and is governed by an equation of state of the form 

p r (r, t) = u>(t)p(r, t) , (46) 

and that the system tends to isotropy for large r. 

For the energy density, we generalize the Mbonye-Kazanas density profile [35| (also utilized by Dymnikova [36|) to 
a reasonable time-dependent model given by 

p(r, t) = p a(<)e- (r/r ° r , (47) 

where po ; r o and n are appropriately chosen constants. Here, the time-dependent function a(t) is chosen so that 
the collapse will asymptote at late times, forming a static star. Note for the sake of clarity that the time-dependent 
function a(t) should not be confused with the junction interface radius introduced in Eqs. ©-([H]). This profile has 
been extremely useful in the investigations of non-singular black holes (i.e., horizons not shielding a singularity) [35| . 
including de Sitter core black holes [36( and, more recently, as a model for gravastars [15( (supplemented with an 
appropriate equation of state). 

The equation of state (|46]) then yields 



Pr (r, t)=u;(t)poa(t)e- (r /^ 



(48) 
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At this stage, it is useful to write the solution to the field equations as follows [37 



-20{r,t) _ 



= 1 - 



8- 



b 2 (t) + / p(f,t)f 2 df 
Jo + 



e 2a(r,t) 

f(r,t) 
Pt{r, t) 



-2/3(r,t) 



exp 



h(t) 



8tt I [ Pr (f,t)+p(f,t)]e 2(3{f - t) fdf 
/o+ 



--j 2b(t)b(t) 

5^ + /) + ( 



p{r, i) r 2 <ir 



2[/3(r,t)-a(r,t)] 



1 



2 a lPr 



(49a) 

(49b) 

(49c) 
(49d) 



respectively, where the explicit coordinate dependence has been dropped in Eq. (|49d[) . and as before, the prime 
denotes a partial derivative with respect to the radial coordinate r, and the overdot a partial derivative with respect 
to the time coordinate t. The functions b(t) and h(t) are two arbitrary functions of integration. For a star, bit) is set 
to zero to avoid a singularity at the center. However, at late time this function need not vanish as it is useful for the 
wormhole configuration. Note that with the prescription of the energy density and radial pressure, the entire system 
of equations may, in principle, be solved for the unknowns. We exploit this fact here. 

To ensure that the star crosses the phantom divide and yet does not collapse for infinite time, the following 
prescriptions are made: 



a(t) = a [(1 



-kot] 



u(t) = Wb - wi (1 - e- klt ) , 



(50) 



where and 1 subscripts denote constant quantities. We now have an infinite family of solutions with the desired 
physical properties. At this stage we should note that the generated space-times tend to Minkowski space-time as 
r — > oo. One would need to therefore cut-off the solution at some r = r* and patch it to an appropriate dark 
energy exterior. However, given how small the energy density and pressures of the currently accelerating universe 
are (when compared to those of an average star) , the asymptotic Minkowski approximation is probably a reasonable 
approximation. 

Surprisingly, for certain values of n, the equations (|49a[) - (|49d[) may actually be integrated to yield an analytic result. 
The expressions are rather unwieldly however so instead we plot the various relevant parameters in figures [5] and [3] 
From the figures it can be noted that the initial inhomogeneity is very small, most pronounced near the center, and 
the space-time almost flat everywhere. There is an inward flow of energy, due to the gravitational collapse as can 
be seen by the negative values of the energy flux plot. At late time, the magnitude of the energy flux decreases and 
asymptotes to zero, indicating the halt of the collapse. At some point during the collapse, the phantom divide is 
crossed and the conditions for possible wormhole formation are established. We show this in figure [4] where it may be 
seen that the NEC violation is most severe in the center. 



IV. TOPOLOGY CHANGE 
A. The Casimir energy approach 

Although (|50[) also allows for an ordinary (positive pressure) fluid at early times, all of the geometries considered in 
the previous sections were modelled so that the evolving parameter u>(r, f)| r =o starts out in the range —1 < uj < —1/3, 
then crosses the phantom divide, and finally ends up in the phantom regime, u> = p r / p < — 1. Once in the phantom 
regime, the negative radial pressure exceeds the energy density, which in principle may imply a topology change. It is 
still uncertain how to obtain this topology change, and if possible, is riddled which difficulties, such as the theoretical 
appearance of closed timelike curves. It is likely that crossing the phantom divide is accompanied by a large quantum 
fluctuation of the metric. Then a crucial question is [jj : what happens when the metric fluctuations become large? 

Concerning the geometry of spacetime undergoing quantum fluctuations, this does not seem to be a source of 
disagreement, but when we turn to the question of whether or not the topology of spacetime undergoes quantum 
fluctuations, the problem becomes more subtle. It was J. A. Wheeler [l9|, [2(| who first conjectured that spacetime 
could be subjected to a topology fluctuation at the Planck scale. This means that spacetime undergoes a deep 
and rapid transformation in its structure. The changing spacetime is best known as spacetime foam, which can be 
taken as a model for the quantum gravitational vacuum. Some authors have investigated the effects of such a foamy 
space on the cosmological constant, for instance, one example is the celebrated Coleman mechanism, where wormhole 
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FIG. 2: Parameters for a sample stellar model presented in section III-C. Plotted on the vertical axes (from left to right, 
starting from the top) are p(r, t), T\.(r, t) p r (r, t), p t (r, t), and finally e 2,3(r ' t) , and -e 2a{r ' t) . 

contributions suppress the cosmological constant, explaining its small observed value (38j |. Nevertheless, how to realize 
such a foam-like space and also whether this represents the real quantum gravitational vacuum are still unknown. We 
can mention some results about topological constraints on the classical evolution of general relativistic spacetimes. 
They are summarized in two points [27| : 

1. In causally well-behaved classical spacetimes the topology of space does not change as a function of time. 

2. In causally ill-behaved classical spacetimes the topology of space can sometimes change. 

From the quantum point of view we can separate the problem of topology change generated by a canonical quan- 
tization approach and a functional integral quantization approach. The Hawking topology change theorem is thus 
enough to show that the topology of space cannot change in canonically quantized gravity [39]. In the Feynman 
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FIG. 4: The figure depicts the null energy condition, NEC = p + p r + 2/ for the specific model with a non-zero flux term 
considered in the text. 



functional integral quantization of gravitation things are different. Indeed, in this formalism, an approach is possible 
to spacetime foam where we know that fluctuations of topology become an important phenomenon at least at the 
Planck scale [40| . However, in our case we can adopt another strategy. In some cases, we can create a one to one cor- 
respondence between topology and the asymptotic energy. In particular, we will consider the Arnowitt-Deser-Misner 
(ADM) energy j4l| as a reference energy. The reason for such a choice is that Eadm > and it is vanishing for 
flat space. Therefore we can think about flat space as the unique reference space to compare a change in spacetime 
associated to the corresponding topology. A trivial example could be the comparison between flat space, where the 
topology is i? 4 and the Schwarzschild space, with topology R 2 x S 2 : they are topologically distinct and possess a 
distinct ADM energy: Efi at — and EschwarzschUd = M . A topology transition from Schwarzschild to flat or 
viceversa, should be necessarily accompanied by a change in ADM energy. In the same manner, we can think that a 
transition from the dark star to the wormhole could be associated to a change in the asymptotic energy, measured by 
the ADM energy, namely if a topology change appears this could be reflected to a change in the ADM energy. The 
way to detect this is simply computed by 

J?DS r?W ormhole / rpDS rpFlat \ ( jpWormhole jpFlat \ > r\ /c-i N 

^ ADM ~ ^ADM — \ ADAf — ADM J _ K^ADM ~ ^ADM) < u - \ 01 ) 

For asymptotically flat spacetimes, the ADM energy is defined as 

Eadm = J (D l h l0 ~ Djh) r j , (52) 

where the indices i, j run over the three spatial dimensions and 

h%j = fjij — (jij , (53) 
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where gij is the background three-metric. Dj is the background covariant derivative and r 3 is the unit normal to the 
large sphere S. However, Hawking and Horowitz fJ|] have shown that the definition (|52")) is equivalent to 

Eadm = 7^ I d 2 x^ (k - k°) , (54) 

07T(_z J ^oc 

where a is the determinant of the unit 2-sphere. k° represents the trace of the extrinsic curvature corresponding to 
embedding in the two-dimensional boundary 2 S in three-dimensional Euclidean space at infinity. In alternative to 
the ADM energy, we can use quasilocal energy to compute such a difference, which is defined by Eq. ([54]) but for 
a finite two sphere. The main reason to use such a definition is that we can extend the surface energy computation 
even to non-asymptoticaaly flat spaces. For this purpose, consider a manifold M composed by two wedges and 
located in the right and left sectors of a Kruskal diagram, respectively and bounded by two three-dimensional 
disconnected timelike boundaries B + and i?_ located in A4 + and M.~ respectively. The quasilocal energy E to t of 
a spacelike hypersurface E = S + U E_ bounded by two spacelike boundaries 5+ and <SL located in A4 + and A4- 
respectively, is given by [HI, [H, [45[ 

Etot = E + — E- . 

More specifically, E to t is defined as the value of the Hamiltonian that generates unit time translations orthogonal to 
the two-dimensional boundaries [H, HI, HEj] • E + and E- are defined as 

(55) 

E - = -1^gSs- d 2 x^(k~k») 
respectively. The trace of the second fundamental form, k, is defined as 

k = -^= (Vhn^ , (56) 

where n M is the normal to the boundaries, and h is the determinant of the metric of E. As an example, consider the 
static Einstein-Rosen bridge, with the metric given by 

ds 2 = -N 2 dt 2 + g yy dy 2 + r 2 (y) dQ 2 , (57) 

where the lapse function N, g yy and r are functions of the radial coordinate y continuously defined on M, with 
dy = dr/-\/l — 2MG/r. The boundaries 5+ and S- are located at coordinate values y — y+ and y = respectively, 

1/2 

and the lapse function is given by \N\ — 1 at both S+ and S-. In this case ro M = (g vy ) 6%. Since this normal is 
defined continuously along E, the value of k depends on the function r, y , which is positive for _B+ and negative for 
B-. See figure [S] for a Penrose-Carter diagram illustrating the boundary locations in a Schwarzschild metric. 
From Eq. (|56[) and Eq. ([57]) . we obtain at either boundary that 

fc=-^, (58) 
r 

where we have assumed that the function r, y is positive for S+ and negative for The trace associated with the 
subtraction term is taken to be fc° = — 2/r for B + and k° = 2/r for B-. As an illustration, consider the case when 
the boundary B + is located at right-hand infinity (y+ = +oo) and the boundary _B_ is located at y_, then 



Etot =M-r 



i 2MG 



(59) 



It is easy to see that E+ and E- tend individually to the ADM mass M when the boundaries 3 B + and 3 S_ tend 
respectively to right and left spatial infinity. It should be noted that the total energy is zero for boundary conditions 
symmetric with respect to the bifurcation surface, i.e., 

E = E+ - E- = M + (-M) = 0. (60) 

Consider now the dark energy star of metric (JlJ and a wormhole defined by the shape function b (r) , with the 
following difference 

Dark energy star (DS) , m (r) with r e [0, +oo) , , 

Wormhole (W) , b(r) with r € [ro,+oo) " ^ ' 
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r = 




r=0 

FIG. 5: A Penrose-Carter diagram illustrating the boundary location in a Schwarzschild metric. M+ and M~ are the two 
wedges, located in the right and left sectors of a Kruskal diagram, respectively and bounded by two boundaries B+ and 
B- located in A4+ and M- respectively. E = E+ US- is a spacelike hypersurface. H+ and H- are the future and past 
horizon, respectively. So (So = H+ n H-) is the bifurcation surface (wormhole throat) and S+ and S- are the two-dimensional 
boundaries of E+ and E_ , respectively. 



Consider also the relation (|5lj) . Thus, by repeating the computation leading to Eqs. (|59|) and (|60[) in the case of 
interest, we get 



,W ,„DS 



k DS = (k-k° 



1 b(r) I 2m (r) 



(62) 



where we are looking at the positive wedge M.+ only. For large boundaries i? 3> r$ and expanding around the throat, 
one obtains 



,w i„DS 



k Ub = — 
r 



b(r )+b' (r )(r-r ) 



-\ 1 



2m (r) 



(63) 



R 



1 - 



b M + fa 7 M (r - rp) + 
2i? 



- 1- 



m(fl) 
i? 



— [r + 6' (r )(r-r )-2m(i2)] 



(64) 



where we have used the wormhole condition at the throat, b(ro) — r$. If b' (ro) =0 and m (i?) is negligible, then we 
recover the ADM mass. Indeed, by integrating on the boundary 2 S+, we obtain 



2G 



M. 



(65) 



If b' (ro) 7^ and m (R) is not vanishing, then the evaluation of the energy depends on a case to case scenario. The 
same discussion can be applied on the negative wedge. As shown in Eq. (|60|) . if we choose boundaries symmetric 
with respect to the bifurcation surface, here represented by the throat ro, we have a total zero ADM- like energy. The 
physical situation looks like a familiar QED physical process 7 — > e + e~: the electric charge is conserved. In our case, 
the charge is the asymptotic energy. Since there is no reason to have an asymmetry in boundaries in the absence of 
external forces, we have to conclude that the classical term is not able to predict the appearance of a wormhole or 
the permanence of a dark star. We are forced to compute quantum effects. The implicit subtraction procedure of 
Eq. (|55|). can be extended in such a way that we can include quantum effects: this is the Casimir energy or in other 
terms, the vacuum energy. One can in general formally define the Casimir energy as follows 



[DM] = E [dM] - E [0] 



(66) 
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where Eq is the zero-point energy, dM is a boundary and Eq [0] represents the zero point energy without a boundary. 
For zero temperature, the idea underlying the Casimir effect is to compare vacuum energies in two physical distinct 
configurations. The extension to quantum effects is straightforward 



E, 



Casirnir 



[dM] = (E [dM] - E 



I classical 



(E [dM] - E Q 



> 1 — loop 



(67) 



In our picture, the classical part represented by the ADM-like energy is vanishing, because of the symmetry of 
boundary conditions. This means that 



E, 



Casimir 



[dM] = (E [dM]-E [0}) 1 _ 



loop 



(68) 



namely Ecasimir is purely quantum. Thus, the Casimir energy can be regarded as a measure of the topology change. 
With this, we mean that, if Ecasimir is positive then the topology change will be suppressed, while if it is negative, it 
will be favored. It is important to remark that in most physical situations, the Casimir energy is negative. Consider 
now the one loop term. We will evaluate it following the scheme of Eq. (|62|) . Thus 



< [dM] - E™ [dM]) = « [dM] E 



1 — loop 



(E o [0]-E™ [dM]), 



-loop 



(69) 



The procedure followed to evaluate Eq. (j6"9"|) . relies heavily on the formalism outlined in Refs. [46J, |47[. The 
computation was realized through a variational approach with Gaussian trial wave functionals. A zeta function 
regularization is used to deal with the divergences, and a renormalization procedure is introduced, where the finite 
one loop is considered as a self-consistent source for traversable wormholes. Rather than reproduce the formalism, we 
shall refer the reader to Refs. 0, [47| for details, when necessary. We can write, 



1 



( E o[9M]-E^[dM]) 1 _ loop = -^ 

+ ( m L ( r ) + m 2,S (Of ln [ — ' 



(m 2 L (r) 



m l,s W) 2 ln 



m 2 L 0) + m\ s (r) 



2.S 



(r) 



Wo 



+ (m 2 L (r) + mis (r)) 2 ln 



(m 2 L (r) 

W L 

m| (r)+ml s (r) 



'1,S 



Wo 



(r)) 2 In 



ml (r) 



n.s 



Wo 




Wo 



(70) 



DS . 



where we have defined two r-dependent effective masses m\ (r) and m\ (r) , which can be cast in the following form 



m\ (r) — ml (r) + 



"i,S 



(r) 



where 



and 



m | ( r ) = m 2 L (r) + m\ s (r) 



b(r) 



ml (r) = -j 1 



™i s (r)=[£zb'(r)-£sb(r)] 



m 2 2 ,s (r)=[^'(r) + ^b(r)] 



(71) 



(72) 



(73) 



respectively. We refer the reader to Refs. [46l |47| for the deduction of these expressions in the Schwarzschild case. 
The zeta function regularization method has been used to determine the energy densities, pi. It is interesting to 
note that this method is identical to the subtraction procedure of the Casimir energy computation, where the zero 
point energy in different backgrounds with the same asymptotic properties is involved. In this context, the additional 
mass parameter /x has been introduced to restore the correct dimension for the regularized quantities. Note that 
this arbitrary mass scale appears in any regularization scheme. Of course b (r) = 2m (r), then we can use only one 
function recalling the different boundary conditions they must satisfy. Generally speaking we can adopt the condition 
m (0) = for the dark energy star and b (r ) — r for the wormhole. Thus, the leading part related to the dark energy 
star close to r — 0, simply becomes 



(EoM-E^ [dM]), 



-loop 



167r 2 r 4 



In 



( 6 



DS 
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On the other hand for the wormhole we get at the throat [5£ 



1 



A(*'M-l) 2 ln 



6' (r ) - 1 



4r 



, !( 6'M + 3) 2 lnf^±^,/ ( 



(74) 



To have an easy comparison with the dark energy star, we make a specific choice for the wormhole shape function. 
We assume that 



then we obtain 



{E?[dM]-E? 8 [8M\) U 



loop g 47r 2 I r 4 \^4rg/i§ 



Moreover, we evaluate the dark energy star term close to r to get 



(EW[dM]-E° s [BM])^ 



1o °p ~ 327r 2 r4 



5 In 



If we choose 



Mo < 



108a 



12. f 6Vi 
— r m ' 



6 In 



6^ 



ro 



(E™[dM]-E° s [8M]), , <0 



(75) 



(76) 



(77) 



(78) 



It is important to remark that the result of inequality (|75|) is valid only for the class of traversable wormholes expressed 
by the shape function ([75]) . To discuss the appearance of different class of traversable wormholes, we need to use 
expression (|74|) inside inequality (|78|) and it is quite evident that this strongly depends on the form of the shape 
function as it should be. It is interesting to note that once this has been created, there is a probability that it will be 
self-sustained [l8j ]. at least for an inhomogeneous u> parameter, like in our case. This means that quantum fluctuations 
related to the Casimir energy play a fundamental part not only for the topology change but even for the traversable 
wormhole persistence. 



B. Morse Index Analysis 



In the classical case there are arguments that if Vo and Vi are compact 3-manifolds, there will exist a space-time 
whose boundary is comprised of the disjoint union of Vq and V\ [23| (see figure [5] for reference.) In relation, Geroch's 
theorem states that if Vq and V\ possess differing topology, then a singularity or closed time-like curves must exist 
somewhere on the manifold 23] . In the realm of wormhole physics one has to generally accept the possibility of closed 
time-like curves unless the kinematics of the wormhole are constrained in some manner [26| . This is true regardless 
of whether there is topology change or not and is simply a consequence of having a wormhole whose mounths may 
move relative to each other. Therefore, in this sense, the issue of closed time-like curves is no more serious a problem 
in the topology changing scenarios than in "standard" wormhole physics. 

Regarding the singularities, even if some singular behaviour exists classically, one may argue that topology changing 
space-times may still contribute to the Lorentzian functional integral approach to quantum gravity where one considers 
the functional integral: 

I = J D[g]e lS . (79) 

Here S is the usual Einstein-Hilbert action. It has been convincingly argued that some of the singularities that arise in 
certain topology changing space-times are extremely mild [23| in the sense that the tetrad becomes degenerate but the 
equations of motion (and the resulting curvature) remain well defined. As well, the Loop quantum gravity approach 
relies on (densitized) tetrads and self-dual connections, and it is know that solutions with classically degenerate tetrads 
yield finite equations of motion also using these Ashtekar variables. Therefore, classically degenerate tetrads are not 
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a) b) 

FIG. 6: Schematics of topology changing space-times via wormhole formation. Figure a) represents topology change via the 
formation of an inter-universe wormhole. Figure b) represents topology change via the formation of an intra-universe wormhole. 
The points p represent the critical point of the topology change. 

necessarily an Achille's heel in this theory. In fact, it is possible that degenerate tetrads may play an important role 
in quantum gravity, [23j | . 

Having established that some of the pathologies associated with topology change are "mild" , the natural question 
to ask is what type of pathology accompanies various topology changes. In this respect, the picture is less clear and 
many studies in the literature are based on case by case bases. Horowitz [23] has convincingly argued that by allowing 
the possibility of degenerate tetrads, topology change is unavoidable. 

Some studies regarding the feasibility of actual topology change rely on the study of Morse functions on the topology 
changing space-times |23[. These studies indicate that manifolds with critical points of Morse index 1 or D — 1 (D 
being the dimension of the manifold) possess causal discontinuities of a severity which are problematic in semi- 
classical analysis (the Borde-Sorkin conjecture [HI). In brief, on the topology changing manifold one constructs a 
Morse function on the metric with Vq and V\ as boundaries [56J. The Morse function, denoted usually by /, possesses 
critical points where d a f = for all values of a. The Morse index, A p at the critical points measures the number of 
negative eigenvalues possessed by the matrix d a dpf at such points. 

In the space-times with Morse index 1 or D — 1, the causal discontinuity causes t he p ropagation of quantum scalar 
fields to become singular somewhere on the manifold, at least in 1 + 1 dimensions [231 ] . This has been used to argue 
that these type of topology changing space-times are suppressed in the sum over manifolds in (|79p . due to the fact that 
they are highly sensitive to small fluctuations. Therefore, such metrics would not make a significant contribution to 
the sum over manifolds due to their combined destructive interference. If this were the case, such topology changing 
spacetimes would be unlikely. 

Figure [6] illustrates a dimensionally reduced schematic of two types of wormhole formation, specifically an inter- 
universe and intra-universe wormhole formation. The critical points of the topology change are denoted by p. In 
such scenarios, the Morse index at p is 1 and such topology change would be suppressed according to the above 
argument in the sum over manifolds approach to quantum gravity. Physically, the causal discontinuity occurs because 
at the onset of wormhole formation, at least two points that were previously not in causal contact suddenly become 
causally connected. (The key point is that this is due to the topology change as opposed to the usual "passage of 
time".) However we hasten to add that at the moment it is not clear that such topology changes are completely 
forbidden, even in the sum over manifolds approach, keeping in mind that small probability is quite different than 
no probability. It is also unknown if the sum over Lorentzian manifolds approach to quantum gravity is indeed a 
valid method to calculate probabilities in a quantum theory of gravity. Although there are now promising candidate 
theories of quantum gravity, it is unknown which, if any, provide the correct methods for calculating properties of 
quantum space-time. 
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V. SUMMARY AND DISCUSSION 

In this work, we have considered time-dependent dark energy star models, with an evolving parameter u) crossing 
the phantom divide, u> = — 1. In particular, we briefly reviewed static and spherically symmetric dark energy stars, 
and further analyzed general solutions of time-dependent spacetimes in detail. Specific time-dependent solutions were 
extensively explored, in particular, the specific cases of a constant energy density, the Tolman-Matese- Whitman mass 
function solution, and a class of models with a non-zero energy flux term, which form from gravitational collapse. 
Once the parameter to evolves into the phantom regime, the null energy condition is violated, which physically implies 
that the negative radial pressure exceeds the energy density. Therefore, an enormous negative pressure at the center 
may, in principle, imply a topology change, consequently opening up a tunnel and converting the dark energy star 
into a wormhole. The theoretical difficulties and criteria for this topology change were discussed in detail, where 
in particular we considered a Casimir energy approach involving quasi-local energy difference calculations that may 
reflect or measure the occurrence of a topology change. Once the topology change has occurred, it is possible that the 
resulting wormhole structures, supported by phantom energy, be self-sustained. As mentioned in the Introduction, 
recent fits to observational data probably favor an evolving equation of state, with the dark energy parameter crossing 
the phantom divide ui — — 1 7]. However, in a cosmological setting the transition into the phantom regime is physically 
implausible for a single scalar field 0] , so that a possible approach would be to consider a mixture of interacting non- 
ideal fluids. One may consider that the time-dependent dark energy star model outlined in this work, is a simplification 
of this possible approach. In fact, recently, static models with two interacting phantom and ghost scalar fields were 
considered, and it was shown that regular solutions exist [28|. It would be interesting to generalize the latter study 
to time-dependent solutions, extending the analysis considered in this work. 

It is interesting to note that the topology change at the center should influence the surface stresses at the thin 
shell, as there is a redistribution of the stress-energy tensor components of the interior solution during the change in 
topology. That this is so may be verified through the conservation identity given by S 1 -^ = [T^e^n"]^ , where [A] + 

denotes the discontinuity across the surface interface S, i.e., [X]* = A+|s — X~\s. The quantity 5j is the surface 
stress-energy tensor at the junction surface S; is the unit normal 4— vector to E; and are the components 

of the holonomic basis vectors tangent to £ (see Refs. [32| for details). Note the dependency of the conservation 
identity on the stress-energy tensor T^, and the right hand side of the conservation identity may also be written as 
5*| ■ = — [or + 2a(a + V)/a\. The momentum flux term, i.e., [Tfj, v e'?^n l ']_ , corresponds to the net discontinuity in the 
momentum flux = U u which impinges on the shell. The conservation identity is a statement that all energy 
and momentum that plunges into the thin shell, gets caught by the latter and converts into conserved energy and 
momentum of the surface stresses of the junction. Now, it may be that the topology change is sufficiently violent to 
disrupture the thin shell. On the other hand, one may also assume that it is sufficiently mild as not to significantly 
affect the stability of the surface layer. 

In analogy to the case outlined in Ref . [22| , where the possibility that inflation might provide a natural mechanism 
for the enlargement of wormholes to macroscopic size was explored, one could imagine that microscopic wormholes 
originated through a topology change, and due to the accelerated expansion of the Universe, these submicroscopic 
constructions could naturally be grown to macroscopic dimensions. For instance, in Ref. (48j the evolution of 
wormholes and ringholes embedded in a background accelerating Universe driven by dark energy, was analyzed. It 
was shown that the wormholc's size increases by a factor proportional to the scale factor of the Universe, and still 
increases significantly if the cosmic expansion is driven by phantom energy. The accretion of dark and phantom energy 
onto Morris-Thorne wormholes [49l |50|. was further explored, and it was shown that this accretion gradually increases 
the wormhole throat which eventually overtakes the accelerated expansion of the universe, consequently engulfing the 
entire Universe, and becomes infinite at a time in the future before the Big Rip. This process was dubbed the "Big 
Trip" [49|, [5(| • However, in the context of the generalized Chaplygin gas, it was shown that the Big Rip may be 
avoided altogether (HHH^]. We refer the reader to Refs. [53[ for more recent details on these issues. In summary, we 
denote these exotic geometries consisting of dark energy stars (in the phantom regime) and phantom wormholes as 
phantom stars. The final product of this topological change, namely, phantom wormholes, have far-reaching physical 
and cosmological implications, as in addition to being used for interstellar shortcuts, an absurdly advanced civilization 
may manipulate these geometries to induce closed timelike curves, consequently violating causality. 

Relative to the topology change issue, a few words are in order. We emphasize that it is still uncertain how to obtain 
this change in topology, and if possible, it is riddled with technical and physical difficulties, such as the appearance 
of closed timelike curves. Nevertheless, it is likely that enormous negative pressures at the center is accompanied by 
a large quantum fluctuation of the metric. The geometry of spacetime undergoing quantum fluctuations does not 
seem to be a source of disagreement in the literature, but the question of whether or not the topology of spacetime 
undergoes quantum fluctuations is more subtle. In the latter, Wheeler conjectured that spacetime could be subjected 
to a topology fluctuation at the Planck scale, where spacetime undergoes a deep and rapid transformation in its 
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structure, resulting in a spacetime quantum foam, which can be taken as a model for the quantum gravitational 
vacuum. Nevertheless, how to realize such a foam-like space and also whether this represents the real quantum 
gravitational vacuum are still unknown. From the quantum point of view we can separate the problem of topology 
change generated by a canonical quantization approach and a functional integral quantization approach. As mentioned 
above, the Hawking topology change theorem is thus enough to show that the topology of space cannot change in 
canonically quantized gravity. In the Feynman functional integral quantization of gravitation things arc different, 
where an approach to spacetime foam is possible where fluctuations of topology become an important phenomenon 
at least at the Planck scale. 
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